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We present the analysis of Fermi Large Area Telescope (LAT) 7-ray observations of HB 21, a mixed- morphology 
supernova remnant. Such supernova remnants are characterized by an interior thermal X-ray plasma, surrounded 
by a wider nonthermal shell emitting at radio frequencies. HB 21 has a large angular size, making it a good can- 
didate for detailed morphological and spectral studies with the LAT. The radio extension is 2°xl?5, compared 
to the LAT 68% containment angle of ~ 1° at 1 GeV. To understand the origin of 7-ray emission, we compare 
LAT observations with other wavelengths that trace non-thermal radio synchrotron, nearby molecular clouds, 
shocked molecular clumps, and the central X-ray plasma. Finally, we model possible hadronic and leptonic 
emission mechanisms. We conclude that 7-rays from HB 21 are likely the result of electron bremsstrahlung or 
proton-proton collisions with dense material due to interaction with the nearby clouds. 
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1. INTRODUCTION 



HB 21 (also known as G89.0+4.7, [Green||2QQ9| ) is a 
mixed morphology (MM) supernova remnant (SNR) 
interacting with molecular clouds. X-rays observa- 
tions with Einstein provide evidence that the super- 
nova (SN) explosion took place in a low-density cav- 
i ty, therefore suggesting a massive stellar progenitor 
(Knoedlseder et al. [1996 ). HB 21 is located in a 
dense environment with a radio sh ape su ggestive of in - 
teractions with the nearby clouds ( |Byun et al. 2QQ6| ). 
Looking at the CO distribution it is possible to dis- 
tinguish two different structures based on the veloc- 



ity of the molecular clouds. |Koo and Heiles 1991 



report the presence of broad CO lines emitted from 
shocked clumps of molecular gas with high column 
densities. The distance to HB 21 was first estimated 



to be 0.8 kpc from ROS AT X-rays observations ( |Leahy 
and Aschenbach [1996 ). Considering the preshock ve- 
locities, the X-ray absorbing column densities, the re- 
lation between photometric distance of Hii region and 
velocity of molecular clouds and the highly polarized 
emission, the distance was later determined to be ~ 1.7 
kpc ([Byun et al.||2QQ6| ). 

High-energy 7-ray sources coincident with 
HB 21 were reported in the 2-years Catalog of 
the Large Area Telescope (LAT) o n board the 
Fermi Gamma-ray Space Telescope (Nolan et al. 
2Q12| ). Three sources are associated with the 



remnant (2FGL J2043.3+5105, J2046.0+4954, and 
J2041. 5+5003) and an additional source lies on 
the edge (2FGL J2051. 8+5054). Due to the large 
apparent size (~2°) this object is well suited for a de- 
tailed morphological study with Fermi LAT. Here we 
present the analysis of ~ 4 years of LAT observations 
of HB 21 and the discussion of the 7-ray emission 
mechanism in light of archival multiwavelength 
observations. In ^2] we describe our morphological 



and spectral characterization of the 7-ray emission 
from HB 21, the sources of systematic errors and 
the strategy adopted to extimate them. In ^3] the 
interpretation of these results is finally discussed. 
This paper is intended to be a summary of the poster 
we presented at Fermi Symposium. For a detailed 
description of the analysis and a more comprehensive 
discussion of the res ults we defer the r eader to the 
paper in preparation Pivato et al. [2013 . 



2. FERMI LAT DATA ANALYSIS 

For this analysis, we use data accumulated from the 
beginning of scientific operations, on 2008 August 4, 
to 2012 June 14, selecting the low-background event 
class P7SOURCE_V6. For the morphological charac- 
terization we use only events with energy >1 GeV to 
profit from the narrower point-spread function (PSF) 
in order to separate the 7-ray emission associated with 
HB 21 from neighboring sources and interstellar emis- 
sion. We then use events down to 100 MeV to deter- 
mine the spectral energy distribution of the remnant. 
Below this threshold the PSF becomes much broader 
than the SNR and the uncertainties related to the in- 
strument response are larger. In both the morpholog- 
ical and spectral characterization we consider photons 
with measured energies up to 300 GeV, but only find 
a significant detection of the source up to energies of 
several GeV due to the limited number of events at 
high energies. 

We perform the analysis in a 10° x 10° region of in- 
terest (Rol) centered at the radio position of HB 21 
(Z=89?0, 6=+4?7). The background is composed of 
diffuse emission and individual nearby 7-ray sources. 
Diffuse emission is taken into account using the stan- 
dard models provided by the Fermi LAT collabora- 
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tion. We include in the background model all the 



point source s present in the 2FGL catalog |Nolan 
et al. 20121 with distances less than 15° from the 



Rol center and which are not associated with the 



SNR. We will discuss in £2.1 the case of the source 



2FGL J2051. 8+5054, which is located at the edge of 
the SNR. The spectral model used for background 
sources is that reported in the 2FGL catalog. Fluxes 
and spectral indices are left as free parameters in the 
fit if the source is within the Rol. Otherwise they are 
fixed to the catalog values. 

Two major sources of systematic errors on the re- 
sults are the uncertainties in the LAT effective area 
and the modeling of interstellar emission. The un- 
certainties in the effective area for the IRFs we use 
are evaluated as 10% at 100 MeV, 5% at 516 MeV, 
and 10% above 10 GeV, linearly varying with the log- 
arithm of energy between those values (Ackermann 
2012] ) . We estimate the error induced in the 



et al. 



characterization of the 7-ray emission spectrum from 
HB 21 by repeating the analysis with two sets of 
modified IRFs where the effective area was upscaled 
or downscaled by its uncertainty. To gauge the sys- 
tematic uncertainties due to the interstellar emission 
model we compare the results obtained using the stan- 
dard model with the results based on eight alternat ive 
interstellar emission models ( |de Palma et al.|[2013 ). 



2.1. Morphological study 

In Figure |TJl we show a count map of the Rol for en- 
ergies > 1 GeV, to visually illustrate the morphology 
of the 7-ray emission in the region. 

We compare the morphology in 7 rays with that 
in the radio (see figure ^) and an X-ray band (see 
figure [2J3) . As shown in figure the 7-ray emission 
is broader than the X-ray (ROSAT) emitting plasma. 
Figure |2^l shows that 7 emission well compares with 
radio shell even if it extends beond the radio shell 
in region where CO clouds are present. Furthermore, 
note that the brightest part of 7-ray emission coincides 
with known shocked molecular clouds (labelled S1,S2, 
and S3 in figure [2^) . 

To characterize the morphology in the 7-ray band 
we approximate the SNR with a disk, for which we 
determine from 7-ray data the best-fit center posi- 
tion and radius. We considered both the cases where 
2FGL J2051. 8+5054 is included as a separate point 
source in the model or removed. The significance of 
the separate source hypothesis is below the thresh- 
old usually required to claim a detection for LAT 
sources so we do not consider 2FGL J205 1.8+5054 
as a separate source and, therefore, we remove it from 
the model. The best fit parameter in this case are: 
/=88?75±0?04, 6=4?67±0?05, r=l?19±0?05. Errors 
reflect statistical uncertainties only. We note that 
the results concerning HB 21 are robust regardless of 



whether 2FGL J205 1.8+5054 is included as a separate 
source in the model or not. 

When the alternative interstellar emission models 
are used, the 7-ray disk is systematically shifted to- 
ward the north-western part with respect to the ra- 
dio shell (with shifts in longitude between 0?19 and 
0?24, and in latitude between 0?06 and 0?09), and the 
disk radius is systematically smaller by 0?18-0?24. We 
note that, as for the standard model, for all the alter- 
native models the disk extends beyond the rim of the 
remnant in coincidence with the western molecular 
cloud. 



2.2. Spectral study 

Modeling HB 21 with the best-fit disk we determine 
the spectrum over the full energy range. We study the 
spectral shape of the 7 emission in the energy range 
between 100 MeV and 300 GeV. We compare a curved 
spectrum 



dN 
dE 



N 



E 



1000 MeV 



-0+/3 ln(£7/1000 MeV)) 



with a power law. We obtain that the curved spetrum 
is preferred at 9a, so now on our spectral shape 
in all energy range will be the log-parabola. Using 
the log-parabolic function, the total energy flux from 
HB 21 at energies above 100 MeV is 9.4 + 0.8 (stat) ± 
1.6 (syst) x 10 -11 erg cm -2 s _1 and the photon flux 
1.48 ± 0.2 (stat) ± 0.4 (syst) x 10" 11 ph cm" 2 s" 1 . 

We also computed the spectral energy distribution 
(SED) in a model-independent way by splitting the 
full energy range in 12 logarithm-spaced bins. We 
show the resulting SED in Figure [3j Systematic er- 
rors related to the modeling of interstellar emission 
are obtained repeating the analysis using the alter- 
native models and extracting the root mean square 
of the variations with respect to the values from the 
standard model. Those values are then summed in 
quadrature with the error due to the LAT effective 
area uncertainties for display. At energies below a 
few hundred MeV the systematic errors are driven by 
the uncertainties related to the modeling of interstel- 
lar emission because the broad PSF reduces the source 
to background ratio. The effective area uncertainties 
have a comparable or larger impact at higher energies. 



3. SPECTRAL ENERGY DISTRIBUTION 
MODELING 

To constrain th e par ticle distribution we simulta- 
neously fit radio (Leahy 2006 ) and 7-ray emission 
from nonthermal electrons and protons. We adopt 



the simplifying assumption that all emission origi- 
nates from a region characterized by a constant mat- 
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Figure 1: a) LAT counts map at energies > 1 GeV. We overlaid the positions of 2FGL sources (crosses for background 
sources and diamonds for the three sources associated with the remnant and 2FGL J2051. 8+5054). The inset in the top 
right corner shows the effective PSF over the energy range considered for a power-law spectral distribution with index 
3.1. b) Background model map (calculated using the fit parameters for the case when HB 21 is modeled using the 
best-fit disk), c) Remaining emission associated with HB 21; overlaid are the positions of the four point sources above 
and the best-fit disk. The pixel size is 0?1 and all maps are smoothed for display with a Gaussian kernel of cr=0?4. 



ter density and magnetic field strength. This sin- 
gle emitting zone is assumed to be equal to the size 
of the remnant derived from the best-fit gamma-ray 
disk. The particle spectra are assumed to follow 
a power-law with an exponential cutoff of the form 
dN/dE ~ T] e ^ p E~ r x ex.p(E/E max ), with the same 
spectral index and energy cutoff for both electrons 
and protons. 

The radio index, a, is related to the 7-ray photon in- 
dex, T, by T = 2a+l for bremsstrahlung or pion-decay 
models, giving T = 1.76 (We note that for this particle 
index, both a cutoff and a spectral break by an index 
of one give an equally valid fit). Modeled spectral en- 
ergy distributions (SEDs) are presented in [4] Model 
parameters are given in Table [TJ We also give the total 
energy of accelerated particles integrated above 1 GeV 
for protons, and above 511 keV for electrons. While 
the chosen parameters are not unique in their ability 
to fit the broadband spectrum, they are representa- 
tive. 



Table I One-Zone Model Parameters 



Model 


r 


Pmax 


n H 


Btot 


Ve/Vp 


w p 








[GeV/c] 


[cm" 3 ] 


[/xG] 




[erg] 


[erg] 


IC 


1.76 


200 


0.1 


3 


1 


1 x 10 49 


4 x 10 49 


Brems. 


1.76 


4 


10 


100 


0.1 


5 x 10 48 


4 x 10 49 


7r°-decay 


1.76 


6 


10 


100 


0.01 


2 x 10 49 


1 x 10 48 



The energetics implied by nonthermal models are 
in line with those determined with other old SNRs 
in a dense environment. Regardless of the domi- 
nant emission mechanism, the energetic requirement 
is for ^fewxlO 49 ergs. The hadronic model indicates 
~ 2 x 10 49 erg in accelerated cosmic ray protons and 
nuclei. Leptonic models require ~ 4 x 10 49 ergs in elec- 
trons. SNR HB 21 is known to be interacting with 
molecular clouds, and a moderately enhanced aver- 
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Figure 2: Emission associated with HB 21 (Fig. [T]:) overlaid with contours from emission at other wavelengths: a) 
Radio emission intensity at 6 cm from Sino-German Telesope |Gao et al.||2011] . Data were smoothed with a Gaussian 
kernel of cr=0?2. The seven contour levels are linearly spaced from Jy beam -1 to 2.0 Jy beam -1 , b) 
Background-subtracted X-ray emission smoothed with a Gaussian kernel of a—0°.25. The three contour levels are 
linearly spaced form 0.36 xlO -3 to 2.13xl0 -3 counts s -1 arcmin -2 . c) Intensities of the 2.6 mm CO line in the 
Local- Arm region. The six contour levels are linearly spaced from 1.5 K km s -1 to 28 K km s -1 . We also show the 
positions of the shocked CO clumps and clouds A and B given in Koo and Heiles 1991 . 



age density would favor hadronic of bremsstrahlung 
models for 7-ray emission. The IC-dominated model 
requires a low density (Nh < 0.1 cm -3 ) to prevent 
bremsstrahlung emission from dominating. Such a 
density is unreasonably low for HB 21. Interestingly, 
HB 21 has an uncharacteristically low luminosity com- 
pared to other 7-ray SNRs known to be interacting 
with molecular clouds. The total luminosity of the 
HB 21 above 100 MeV at a distance of 1.7 kpc is 
4 x 10 34 erg s -1 , while most interacting SNRs, such 
as IC 443 have luminosities of 10 35 erg s -1 or greater. 
It is possible that the lower luminosity of HB 21 is 
due to a lower total mass of the molecular clouds with 
which it is known to be interacting. 



4. CONCLUSIONS 

We analyzed the 7-ray measurements by the Fermi 
LAT in the region of HB 21, a mixed- morphology 
SNR. We detect significant emission (~29a associ- 
ated with the remnant. The emission is best modeled 
by a disk centered at (l,b) = (88?62±0?05,+4?79±0?06) 
with a radius r=l?19±0?06 (statistical uncertainties 
only), so it is well resolved by the LAT for energies 
greater than 1 GeV. The 7-ray emission extends over 
the whole area of the non-thermal radio shell, larger 
than the X-ray emitting thermal core. The emission in 
7-rays may extend beyond the radio shell in a region 
rich of interstellar matter in the north western part 
of the SNR. Furthermore, the brightest 7-ray emit- 
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Figure 3: Spectral energy distribution (SED) of HB 21. 
The line shows the best-fit log-parabola model, the light 
filled area shows the statistical error band and the gray 
area shows the systematic uncertainties. For the points 
error bars correspond to statistical errors only, while lines 
show the larger systematic errors. 95% confidence- level 
upper limits are given for energy bins where the TS of 
the source is < 9. 




Figure 4: SED models for which IC (top), 
bremsstrahlung (middle) and 7r°-decay (bottom) are the 
dominant emission mechanism. In each model the radio 
data are fit with a synchrotron component, shown as a 
black curve. The individual contributions of tt° -decay 
(long dashed), bremsstrahlung (short dashed), and IC 
emission from CMB (dotted) are shown. The sum of the 
7-ray emission is shown by the solid curve. The leptonic 
and hadronic components are colored blue and red, 
respectively. 



ting region coincides with known shocked molecular 
clumps. Both results are suggestive that collisions 
of shock-accelerated particles with interstellar matter 
are responsible for the observed 7-ray emission. 

The spectrum is best modeled by a curved func- 
tion, indicative of a cutoff or break in the spectrum of 
the accelerated particles, typical of middle-aged SNRs 
in a dense interstellar environment. The total 7-ray 
luminosity of HB 21 above 100 MeV is estimated to 
be ~ 4 x 10 34 erg s -1 , fainter than other SNRs in- 
teracting with MC detected by the LAT. This can be 
explained by the lower mass of the molecular clouds 
supposed to be in interaction with the remnant. The 
7-ray emission may be dominated either by 7r° de- 
cay due to nuclei or by Bremsstrahlung from energetic 
electrons. An IC origin is disfavored because it would 
require unrealistically low interstellar densities in or- 
der to prevent Bremsstrahlung to dominate. Based 
on the most likely values for the ISM densities over 
the volume of the remnant, in the hadronic-dominated 
scenario accelerated nuclei contribute a total energy of 
~2xl0 49 ergs. In the leptonic-dominated scenario we 
need accelerated electrons for ~4xl0 49 ergs, and also 
nuclei for ~10 49 ergs. Such energy densities match 
those required to make SNRs the do minant source of 
the Galactic CRs [e.g . |Q'C. Drury||2012| . Recently, 



Reichardt et al. |2012| performa an analysis on HB 21 
using 3.5 years of Fermi LAT data. The conclusions 
are in agreement with those presented in this proceed- 
ing if we consider that w e adopt a distance from HB 21 
of 1.7 kpc (according to Byun et al. |2006| ) instead of 
0.8 kpc adopted in earlier works and based on the as- 
socition with the Cygnus OB7 complex. In particular 
our values of flux and luminosity of the SNR are in 
agreement. Finally they assume a density of 60 cm -3 
which is the maximum density estimation. Insted, we 
consider both the minimum and the maximum possi- 
ble desities obtaining that the energy accelerated par- 
ticles is ~ 0.3 — 4 x 10 49 erg, regardless of whether 
hadronic or leptonic emission is dominant. Finally, 
in our work also systematic effects due to interstel- 
lar emission model and effective area are taken into 
account. 
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